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ABSTRACT 

We study dissipation process of magnetic fields in the metallicity range 0 — IZ© for contracting 
prestellar cloud cores. By solving non-equilibrium chemistry for important charged species including 
charged grains, we evaluate the drift velocity of the magnetic-field lines with respect to the gas. We 
find that the magnetic flux dissipates in the density range 10^^cm“^ ^ ^ for the solar- 

metallicity case at the scale of the core, which is assumed to be the Jeans scale. The dissipation density 
range becomes narrower for lower metallicity. The magnetic held is always frozen to the gas below 
metallicity < 10“^ — 10' depending on the ionization rate by cosmic rays and/or radioactivity. 

With the same metallicity, the dissipation density range becomes wider for lower ionization rate. The 
presence of such a dissipative regime is expected to cause various dynamical phenomena in protostellar 
evolution such as the suppression of jet/outflow launching and fragmentation of the circumstellar disks 
depending on the metallicity. 

Subject headings: early Universe—metal-poor stars 


1. INTRODUCTION 

The first episode of star formation makes a large 
impact on the subsequent thermal and chemical evo¬ 
lution of the universe by initiating such events as 
metal enrichment and reionization of the intergalac- 
tic medium. The first stars, so-called population 
III stars, are formed from the primordial pristine 
gas, consisting only of hydrogen, helium and a trace 
amount of deuterium and lithium. Theoretical stud¬ 
ies predict that they are typically very massive ^ 


lOOMo) (e.g. iBromm et al.ll2002l: lOmn 

[kai & Palla 

to 

o 

o 

McKee & Tan 1 120081: iHosokawa et al. 

120111: ISiisa 

2013 


Hirano et al.ll2Q14t iSusa et aI]l2Q14l ). and some of them 


end their li ves as supernova explosions after a few mil¬ 
lion y ears (jUmeda fc Nomoto I I2QQ2I: iHeger fc Wooslevi 

Hooi. 

Following the explosions, metal-enriched ejecta is 
mixed with the ambient material and may recollapse 
to form the next generation of stars (e.g., Ritter et 
al. 2012). Metallicity of the newly formed stars re¬ 
flects th e degree of mixing of the ejecta w ith the pris¬ 
tine gas (jWise et al.ll2Q12l: iChen et al.ll2Q14l ). A gas with 
metals cools more efficiently owing to the higher radi¬ 
ation emissivity as well as higher sustainability against 
stellar radiative feedback. Since the cold environment 
thus achieved is favorable for gravitational collapse and 
fragmentation of star-forming clouds, vigorous star for¬ 
mation might be induced by the metal enrichment . 
In fact, in som e numerical models (jAhn et al.l I2Q12I: 
IWise et ^ l2Q14l ). second-generation stars formed from 
metal-enriched gas, rather than the first stars, are the 
dominant radiation sources for the cosmic reionization. 
Knowledge of star formation in low-metallicity environ¬ 
ments is a clue to understand the cosmic structure for¬ 
mation. 

Thermal and dynamical properties of low-metallicity 
star-forming clouds has been investigated by a num- 
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ber of semi-analytic (lOmukail 20001: ISchneider et al.l 

120021: lOmnkai et, all 1200.1 12010: I.Ti et all 1201411 

and 

numerical studies (iBromm et al.l l2001l: iMachidal 

2008: 

Clark. Glover. & KlessenI 20081: iDoocke et al.l 

2011, 

20131). With slight metal enrichment of 10“^ — 10“ 

-4 y 

^05 


the cooling by dust grains causes rapid cooling at high 
enough density (> lO^^crn”^) where the Jeans mass is 
small (< IMq), thereby triggering fragmentation of the 
clouds into small pieces (I Omukail 1200(1 iSchneider et al.l 
2QQ2I. l2QQ3l: lOmnkai et alT 2QQ5I: iSchneider et al.l I2QQ6I: 
Omnkai et al.l I2Q1QI: ISchneider et al.l l2Q12f ). At this 
critical metallicity, whose exact value depends on the 
nature and amount of dust in the gas, the transition 
of characteristic stellar mass fro m massive to low-m ass 
ones is expected to occur (but see lFrebel et al.l[2QQ7[ ). 

Most studies so far on low-metallicity star forma¬ 
tion have not taken account of magnetic fields assum¬ 
ing that only weak fields are present in the early uni¬ 
verse. Those seed fields are then amplified by the as- 
trophysical dynamo, including by galactic differential ro¬ 
tation, to the present /iC level in the local ISM, com¬ 
parable in energy density with thermal and turbulent 
energies. In present-day star formation, the magnetic 
fields indeed play colorful roles such as launching pro¬ 
tostellar jets and outflows and tr ansffering angulp mo - 
mentum in protostellar disks (e.g. iMachida et aTll2QQ8f) . 
Recently, it has been recognized that, even in the pri¬ 
mordial environment, magnetic fields may be ampli¬ 
fied by so-called small-scale turbulent dynamo and the 
field strength may reac h a dynamically signi fi cant level 
(iSchleicher et al.l I2Q1QI: l^chober et al.l l2Q 121: iSur et al. 

B iFederrath et al.l[2QTTI: iTurk et ^120121: iLatif et ^ 

. Once the first stars are formed, magneti c fields 
can b e generated and amplified in the stars (e.g. iSoruitl 
l2QQ2f ). the circumstellar disks (jShiromoto et al.l 120141)^ 
and e ven the neighbori ng prestellar cores (jAndo et al.l 
|2Q1QI : iDoi fc SusgJ l2Qllf ). Finally, at their supernova 
explosions, those fields are dispersed into interstellar 
and intergalactic medium. Therefore, when the second- 
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generation stars are formed from re-collapsed gas, mag¬ 
netic fields are believed to be already present. Although 
proper assessment of the field strength in such environ¬ 
ments is still elusive, it is worthwhile to study the effects 
of magnetic fields on low-metallicity star formation. 

Recently, magneto-hydro dynamical (MHD) simulation 
of star formation in low-me tallicity, magnetized clouds 
has been carried out by iPeters et al.l (j2Q14f ) under 
the ideal assumption, i.e., perfect flux freezing to the 
gas. They found that the magnetic field generated via 
the small-scale dynamo stabilizes the protostellar disk 
against fragmentation, counteracting destabilizing effect 
due to the enhanced cooling by metals. However, it is 
still unclear in what metallicity range the ideal MHD 
approximation is justified. For a metal-free gas, it is 
known that the coupling between the field and gas is so 
tight that ideal MHD is a good assumption for the Jeans 
scale (|Maki fc Susal[2QQ^ 12QQ7D . On the other hand, in 
the present-day solar metallicity gas, this approximation 
breaks down and the fields dissipate in the density range 
10^^ — 10^^cm“^. As a result, the actual protostellar 
disk formation process can be different from that with 
the ideal approximation, which overestimates magnetic 
effects on dynamics. In view of this, we here investigate 
the coupling between the gas and magnetic field, or in 
other words, resistivity of the gas, at various metallic- 
ities between 0 and IZq. Since it is the abundance of 
charged particles that controls the resistivity, we solve 
non-equilibrium chemical reaction network for prestellar 
cores, which are assumed to be collapsing dynamically. 
We then examine the validity of the flux-freezing condi¬ 
tion during the prestellar collapse. 

This paper is organized as follows: in Section 2, we 
describe the thermal/chemical evolution of the prestel¬ 
lar cores, which is used as a background model for the 
analysis of the magnetic-flux dissipation described in Sec¬ 
tion 3. Finally, in Section 4, we summarize our results 
and discuss its implication on the star formation in low- 
metallicity environments. 


2. THERMAL AND CHEMICAL EVOLUTION OF 
COLLAPSINC PRESTELLAR CORES 

2.1. Numerical modelling 

We here employ one-zone model developed by 
lOmnkai et al.l (j2QQ5h , in which a spherical cloud with its 
core radius of the Jeans scale Rj is assumed to collapse 
in a runway, self-similar fashion, with some modification 
in chemical reactions. The calculated quantities in this 
model are those at the center. The density p of the cloud 
increases as 


dp 

dt 




Tw, 


where the free-fall time 


= 



( 1 ) 


( 2 ) 


G the gravitational constant, and / the ratio of the pres¬ 
sure gradient force to gravity at the cloud center, which 
is fitted approximately as a function of the effective ra¬ 
tio of specific heat 7 = d In p/d In p (p denotes the gas 


pressure): 


r 0 7 < 0.83 

f = \ 0.6 + 2.5(7 - 1) - 6.0(7 - 1)^ 0.83 < 7 < 1 

I 1.0 + 0.2(7 - 4/3) - 2.9(7 - 4/3)^ 7 > 1. 

( 3 ) 

The thermal evolution is followed by solving the energy 
equation for internal energy per unit mass e: 


dt dt \pJ 


-A, 


net 5 


( 4 ) 


where Anet is the net cooling rate of the gas. The equa¬ 
tions above are supplemented with the equation of state 


P = (7ad - l)pe, 


( 5 ) 


where 7ad is the adiabatic index. The temperature T and 
energy density e are related with 


1 fcsT 

7ad - 1 prrip 


(6) 


where /cb is the Boltzmann constant, p the mean molec¬ 
ular weight, and rup the proton mass. For the net cool¬ 
ing rate Ane t, in addition t o the o riginal cooling/heating 
processes in lOmnkai et al.l (j2QQ5[ ). which include the ra¬ 
diative cooling by [CH], [Cl], [01], H2, HD, CO, OH, 
H2O lines, and by continuum from the primordial gas 
and dust, and cooling and heating associated with chem¬ 
ical reactions ( see appendix), we also consider the heat¬ 
ing due to ionization by cosmic rays an d radioactivity. 
The io nization heating rate is taken from lWolfire et alTI 
(|1995[) with assumption that the heat deposition rate as¬ 
sociated with ionization by the decay of radioactive el¬ 
ements (REs) is the same as by cosmic rays. We do 
not include the heating due to the dissipation of mag¬ 
netic fields. It has little impact on the thermal evolu¬ 
tion since the thermal energy is always larger than the 
mag netic energy on avera ge under realistic circumstances 
(e.g. iTomida et al.|[2QT^ . 

The chemical fractions of the coolants are followed by 
solving the non-equilibrium chemical network. We add 
14 new species Li, Li+, Li^+, Li^+, Li“, LiH, LiH+, M, 
M+, G, G+, G-, G2+, G^-, which can be irnportan t 
charge providers, to the model of lOmukai et al.l (|2QQ5F 
Here, M and G stand for the metallic elements and 
grains, respectively. The re actions related to L i and its 
ions/molecules are listed in iBovino et al.l (l2nnh and in¬ 
cluded in the present calculations. The initial number 
abundance of Li relative to hydrogen, yu = ^Li/^H, 
is taken as 4.8 x 10“^^, the observed ISM value in 
Small Magellanic Glo ud, consistent also with the big 
bang nucleosynthesis (|Howk et al.l I2Q12I: iGalli fc Fallal 
|2Q13[ ). M and M+ represent atoms and ions, respec¬ 
tively, of the metallic elements Na, Mg, Al, Ga, Fe and 
Ni, which are possible major electron providers. Since 
they all have similar rate coefficients for the charge 
transfer reaction with other ions as well as for radia¬ 
tive recombination with electrons, we ignore their dif¬ 
ference and treated collectively by summing up their 
elemental abundan ces. The number fraction o f M is 
Pm = 1.68 X 10“^ (jUmebavashi fc Nakanol[T99Qf ) at the 
solar metallicity and is proportionally reduced in lower 
metallicity cases. The co efficients of M and M+ r elated 
reactions are taken from lUmebavashi fc Nakanol (|199Ql ) 
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and iPrasad fc HuntressI (|198Qf ). Reactions relevant to 
the dust grains are crucial to determine the ionization 
state. We consider five ionization states (G, G^, G“, 
G^+, G^“) of the dust grains. We assume the same mas s 
fraction and size distribution as in lOmukai et al 1 (1200^ . 
Mass fraction of the dust is 0.9 39 x 1Q~^ below the w ater- 
ice evaporation temperature (jPollack et al.l 119941). and 
the s ame size distribution is assumed as (jMathis et al.l 

[ 1973 ): 

dugr f 5 X 10“^/im < a < 1/im 

da ^ \ 1/im < a < 5/im. ^ ^ 

The evaporation of each component of the dust takes 
place instantaneously at its vaporizatio n temperature 
(|Omnkai et al.l I2QQ5I: iPollack et al.l 119941) . The collision 
rates between grain-charged part icles and grain-g r ain ar e 
calculated by eqs.(3.1)-(3.5) of iDraine fc Sntinf (jl987[ ). 
averaged over the size distribution (eq. [7)). 

Ionization by cosmic rays and radioactivity controls 
the ionization degree in the clouds. The cosmic rays 
ionize atoms and molecules either directly or indirectly 
by high-energy photons emitted at the direct ionization 
(described as GRPHOT). The dire ct ionizat i on rat es by 
cosmic rays are same as those of lOmnkail (j2Q12f ) , ex- 
cept for the rates of M and H CO, which are taken from 
lUmebavashi fc Nakanol (IT^ . The GR PHOT rates of 
M an d HGO are taken from UMIST2 (jMcElrov et al.l 
l2Q13f ). The attenuation of cosmic rays is considered for 
ionization rate by cosmic rays (cR' 

CcR = vCcR,o exp (^) 

where CcR,o is the cosmic-ray ionization rate in the local 
ISM and A the attenuation length. We use Ccr,o = 1 x 
10“^^s“^ and A = 96 g cm“^ in the present calculations. 
In the above, a parameter 77 is introduced to specify the 
amount of ionization sources. The core radius is given 
by the Jeans length 

<»> 

REs also ionize the gas by emitting gamma rays at 
their decay. We classify REs into two categories by the 
decay time and treat their contribution to ionization rate 
differently. We assume the amount of long-lived REs, e.g. 

with decay time tdecay ^ IGyr proportional to the 
total metallicity since the long-lived REs accumulate in 
the ISM with little decay as the chemical enrichment pro¬ 
ceeds. By using the value in the solar neighborhood fo r 
the solar metallicity case (|Umebavashi fc Nakanol|2QQ^ . 
its contribution to the ionization rate is given by 

CRE,long = 1-4 X 1O“^^S“^Z/Z0. (10) 

The short-lived REs, e.g.^^Al, have larger effects than 
the long-lived REs since the gamma-ray intensity is in¬ 
versely proportional to the decay time. The short-lived 
REs decay in shorter timescale than the age of the uni¬ 
verse at high redshift and hence do not accumulate in 
the ISM. Without detailed knowledge of their amount, 
we here simply assume that it is also proportional to the 


same ionization parameter 7/ introduced above: 

CrE, short = 7.6 X 10“^^s“^7/. (11) 

In summary, the total ionization rate by the cosmic 
rays and the REs can be written as 

C = CcR + CrE, short + CRE,long- (12) 

The first two terms are proportional to the parameter 
77 and the last term is to the metallicity. In this paper, 
we consider four cases with 77 = 0, 0.01,1,10, which are 
denoted as models 1, 2, 3 and 4, respectively. Model 1 
(77 = 0) is the case without the ionization either by cos¬ 
mic rays or short-lived REs. In particular, for Z = 0, no 
ionization source is present in this model, which can be 
regarded as the primordial pristine environment. Model 
2 (77 = 0.01) has 100 times smalle r ioniz ation rate than 
the local value. iStacv fc BrommI (|2007[) estimated the 
cosmic-ray intensity is about this value in first galaxies 
by summi ng up the contribution from population III su¬ 
pernovae. iNakanchi et al.l (|2014D also assessed Ccr,o at 
6 ^ ^ ^ 1^5 spreading over two orders of magnitude 
around 10“Model 3 (77 = 1) has the cosmic-ray 
intensity the same as in our Galaxy and allows us direct 
comparison wit h previous studies on present-d ay star for¬ 
mation such as iNakano fc Umebavashil (|1986[) . In model 
4 (77 = 10), with supernova explosions ten times more 
frequent than in our Galaxy, the environment would be 
similar to starburst gal axies at z < 5 with much intense 
star-forming activities (jLackill20M ). 

We calculate nine cases with different metallicities 
Z/Zq = 0,10“^10“^10“^ 10“^, 10“^ 0.01, 0.1,1 for 
each ionization model. Numerical convergence of non¬ 
equilibrium chemical reactions occasionally becomes very 
slow at high densities. To avoid such difficulty, we switch 
to the equilibrium Saha equations of the H and He sys¬ 
tem for tth > 10^^cm“^ neglecting all the metal species. 
In some cases, the convergence becomes extremely slow 
at lower densities and we are forced to switch to the 
equilibrium calculation at 10^^cm“^ < < 10^^cm“^. 

In metal-enriched cases, this treatment causes a discon¬ 
tinuity in the fraction of some charged species around 
10^^cm“^. Those discontinuities comes from the fact 
that we are still omitting some reactions required to 
achieve the chemical (i.e., Saha) equilibrium such as the 
ionization by thermal photons, etc. Note that to reach 
the correct Saha equilibrium at high enough density, both 
the forward and reverse reactions must be included for 
all the relevant processes. In addition, we do not include 
the species other than e“, H, H+, He, He+, He++ in solv¬ 
ing the Saha equations, which can be part of the reason 
of the discontinuity. Q 

However, this does not affect significantly our conclu¬ 
sion on the density range of magnetic field dissipation. 
This is because the discontinuity appears at high enough 
density where the charged dust is already vaporized. In 
the absence of charged grains, electrons dominate the 
electric current carrier, which is more or less properly 
estimated by the Saha equation in such a high density 
regime. 

^ Vaporization of the grains also generate discontinuities at 
~ because we model the vaporization of grain materials 

by the threshold temperatures. Hence the discontinuities around 
~ 10^®cm“^ is reasonable. 
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2 . 2 . Thermal evolution 

In this section, we describe the thermal evolution of 
the clouds for the four ionization models. The panels 
in Fig(T] show the evolutionary tracks on the uyl — T 
plane for models 1 - 4 (77 = 0,0.01,1 and 10). Different 
curves in each panel are for different metallicities. Below 
^ 10 ^^cm“^, the temperature decreases with increas¬ 
ing metallicity owing to the larger amount of coolants 
for a given ionization parameter. After the central part 
of the cloud becomes optically thick to the dust ther¬ 
mal radiation, all the evolutionary tracks with different 
metallicities co nverge to a single track for > 10 ^^cm“^ 
(| 0 mnkaill 2000 [) . 

With larger ionization parameter, the temperature 
tends to be lower for the same metallicity and density 
range. For example, with a small amount of cosmic-ray 
ionization (model 2 , 77 = 0 . 01 ), the temperature is lower 
than in the no-ionization case (model 1 , 77 = 0 ) in the 
range 10^ —10^cm“^. This tendency is more clear in mod¬ 
els with higher ionization rate, i.e., models 3 (77 = 1) and 
4 (77 = 10). This is because ionization promotes H 2 for¬ 
mation via the electron catalyzed reaction. By enhanced 
H 2 cooling, the gas reaches to lower temperature, which 
is favorable for HD formation. Hence, the temperature 
plummets further by the HD coo ling with its excitation 
energy four times lower than H 2 (jStacv fc BrommlDOOTl: 
iNakanchi et al .11201^ . With Z > IO^^Zq, other cooling 
agents such as C and O are more important than HD and 
the temperature evolution is little affected by the differ¬ 
ence in ionization rate. In all cases, the temperature in 
high-density regime {nu > 10 ^^cm“^) does not depend 
on the ionization parameter because of cosmic-ray atten¬ 
uation as well as dominance of dust cooling. 

2.3. Chemieal abundanees 

Since we are interested in the resistivity of the gas, we 
plot the abundances of 7 species (e“, G^, H+, HCO+, 
M+, Li+), which can be important charge carriers, as a 
function of number density in Fig. [2l In this Figure, 
16 panels correspond to models with 4 different metal- 
licities (Z/Zq = 1,10“^, 10“^ and 0) and 4 ionization 
rates (77 = 0,0.01,1, and 10). In all cases, e“ and H+ 
are the dominant charged particles in the low-density 
regime, but their recombination proceeds with increas¬ 
ing density. The charged grains eventually become 
the dominant charged component at some density except 
in the zero-metallicity cases (i.e., the rightmost column). 
At even higher densities (> 10 ^^cm“^), the vaporiza¬ 
tion of grains quenches the recombination on the dust 
surface, and thus the ions and electrons dominate the 
charge again. The behavior in the zero-metallicity cases 
is quite different. The initial decrease of the electron frac¬ 
tion by the recombination with H+ is followed by a floor 
set by Li+. Hence, the ionization degree never becomes 
lower than 10“^^ for all 77 values. Note that this behav¬ 
ior of the ionization degree in zero-metallicity clouds has 
been reported for smaller ranges of ionization pa rameter 
(jMaki fc Snsall2QQ3. l2QQ7l : iSchleicher et al.l[ 2 QlQf ). 

Next we see the difference among the ionization mod¬ 
els. In the 77 = 0 case (model 1; first row), the recom¬ 
bination proceeds much faster than in non-zero 77 cases 
(the second to the fourth rows). The electron fraction 
drops exponentially as a function of the density for the 


77 = 0 cases (model 1 ), while it decreases more slowly in 
other cases as a result of the cosmic-ray/RE ionization. 
At a fixed density, the electron fraction increases with 
the ionization rate. The higher electron fraction means 
the less importance of charged grains: the density range 
where the grains dominate the charge shrinks toward in¬ 
creasing 77 . For cases with 77 = 1 and 10 (models 3 and 4) 
and with IO^^Zq, the charged grains never dominate the 
charge because of the small amount of dust and relatively 
high ionization rate. 


3. DISSIPATION OF MAGNETIC FLUX 

Dissipation of magnetic fields in a gas can be studied 
by estimating the ambipolar diffusion rate and Ohmic 
resistivity. Here we ernploy t he formalism developed by 
iNakano fc Umebavashil (jl986f ) , which allows us to handle 
the two dissipation processes simultaneously. 

The drift velo city of magnetic field lines rel ative to the 
gas is given by (jNakano fc Umebavashil [l986f ): 


where 


VBx 




U X B), 


A = Af + Al 




pvTyUJZ 


Ao = 


' V V 


(13) 

(14) 

(15) 

(16) 


Here c is the speed of light, j and B the electric current 
and magnetic flux density, pj^ and ujjy are the density and 
the cyclotron frequency of the charged particle z/, respec¬ 
tively. The subscript x denotes the direction parallel to 
j X B. The viscous damping timescale of the relative 
velocity between charged (z/) and neutral (n) particles is 
given by 

Tu = - — - (17) 

Pun'^u'^n ^O'V^ i/Yi 


where '^n are the reduced mass of particles u and 

n, number densities of u and n species, respectively. The 
momentum-transfer rate coefficient via th e collision be¬ 
tween z/ and n, <av>jyn^ is evaluated as in iPinto fc Gallil 
(|2008f) and averaged over the dust size distribution 0 . 
Here, we replace the term (jf x B)^ in equation ([13]) by 
the mean value of magnetic force where B is 

the mean magnetic field s trength and R the scale lengt h 
that we are interested in (|Nakano fc Umebavashil (19861 ). 
which is taken as the Jeans radius Rj (eq. 

Relative importance of the Ohmic loss to the ambipo¬ 
lar diffusion is determined by the quantity \rj^ujy\. If 
it is less than unity, the Ohmic loss is more impor¬ 
tant than the ambipolar diffusion, and vice versa. Ap- 
proximate expressions for tw o limits can be given as 
(|Nakano fc Umebavashil ll 9861 ): 


VBx 


AiracR 
Tj B^ 

^ Pi AttR 


{\rjyUJjy \ <C 1 Ohmic loss) 

{\ryUJy \ ^ 1 ambipolar diffusion) 

(18) 
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Fig. 1.— Temperature evolution of the prestellar cores for different metallicities as a function of the number density. Four panels (a)-(d) 
correspond to cases with four different ionization parameters 77 = 0, 0.01,1 and 10, respectively. The values of metallicity are indicated by 
numbers in the Figure. In panel (a), the adiabatic temperature gradient for 7 ad = 7/5 is shown by the dashed line. 


where dc is the conductivity defined as 




^ rriu ’ 


(19) 


and Qi, and rrijy are the charge and the mass of a charged 
particle z/, respectively. The subscript i represents the 
dominant ion species. 

We compare the drift velocity vbx (eq. [18]) with the 
free-fall velocity ixff defined as 

/2GM 

uft = y ^ ■ ( 20 ) 

For the mass and radius of the prestellar core M and R, 
we use the instantaneous Jeans mass (Mj = AirpRj/S) 
and radius (eql9j), respectively. If vbx > i^ff, the magnetic 
field lines move away from the core by their own tension 
before the significant collapse, i.e., the magnetic fields 
dissipate. On the other hand, for vbx < i^ff, the magnetic 
fields hardly dissipate and the field lines can be regarded 
frozen to the cloud. 

FigsISM] show the ratio VBx/uf^ on the ub -B plane 
for four cases with different ionization parameters r] = 
0,0.01,1 and 10 (i.e., models 1-4, respectively). In each 
figure, four panels show the cases with metallicities 1^©, 
10“^Z©, 10“^Z© and 0. In the red regions in those pan¬ 
els, the condition vbx/u^^ ^ 1 is satisfied, i.e. the mag¬ 
netic flux dissipates from the gas during the collapse. 

We describe the results for the r] = 1 case (Fig El model 
3) as the fiducial model in the followings. The qualita¬ 
tive behavior is common in all the models although with 
some quantitative differences. First, in relatively metal- 
enriched cases Zq and 10“^^©, i.e., top two panels (a) 
and (b), a red strip parallel to the ^-axis appears in the 
density range 10^^ — 10^^cm“^, i.e., the magnetic fields 
dissipates between this interval irrespective of the held 


strength B. This is due to the Ohmic loss, for which the 
drift velocity vb^ does not depend on B (see eq. EHl), and 
so the dissipation occurs regardless of the held strength. 
This strong dissipation comes from the dominance of the 
grains in charge in this density range as seen in Fig. O 
The grains have much larger inertia than the electrons, 
causing large resistivity to the current. Existence of such 
a dissipation density range has al ready been pointed out 
by iNakano fc Umebavashil (j 19861 ) for the local ISM con¬ 
dition. Our study reveals that it continues to exist even 
for lower metallicity gases. In extremely metal-poor cases 
(see panels c and d ) for Z = 10' and 0, respec¬ 
tively), the dissipation range becomes very narrow for 
Z = 10“^Z© (panel c) and almost disappears Anally for 
Z = 0 (panel d) because the grains never dominate the 
charge owing to their small amount (see Figl2j) and so 
the Ohmic dissipation does not work. In addition, the 
recombination of charges on the grain surface does not 
proceed efficiently at such low metallicity, and the abun¬ 
dance of electrons and ions remains relatively high. In 
particular, in the Z = 0 case (panel d), Li+ remains in 
the gas phase because of the complete absence of dust 
grains, which would absorb Li+ otherwise. As a result, 
the electron supply from Li+ gives a floor to the elec¬ 
tron fraction, which keeps the resistivity very low and 
diminishes the dissipation range. 

In all the panels in Figs El - El dissipation regions by 
the ambipolar diffusion are present at the upper left (low 
density and strong held) corners. From equations ([T8j) 
and (jSQ]) at the Jeans scale, the boundary of the dissipa¬ 
tion region defined as 'Z^bx/'^f = 1 is given by 

B oc (21) 

roughly consistent with the slope of the contour. Note 
that the magnetic energy density is comparable to the 
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Fig. 2 .— Abundance of major charged species as a function of number density. Results for ionization parameter 77 = 0,0.01,1 and 10 
(models 1-4) are shown from the top to the bottom rows, and for metallicities Z = Z©, 10~^Zq, Zq,0 from the left to the right 
columns. The negatively charged species e“ and G“ are shown by thick solid curves, while the positive ions (H+, G“^, HGO“^, M+, Li+) 
are plotted by thin solid lines. Note that the discontinuities seen at ~ 10^®cm“^ for metal-enriched cases are artifacts caused by our 
switching from non-equilibrium to equilibrium chemistry. 
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Model 1 (ri=0) logvBx/uff 
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lognH[cm'^] lognH[cm'^] 


Fig. 3.— The ratio vbxI^s on the nn — B plane for the 77 = 0 cases (model 1). The numbers indicated on the contours denote its 
logarithmic values. Four panels (a)-(d) correspond to different metallicities (a)Z/Z 0 = 1, (b)10“^, (c)10“®, and (d)0, respectively. The 
white curves denote Bcr{see text). 























































Fig. 4.— Same as Fig[3l but for the ionization parameter 77 = 0.01 (model 2 ). 


gravitational energy density, and the collapse is prohib¬ 
ited in this parameter range. The critical magnetic flux 
Bcr, which is defined by the relation 

= pGMj/Rj, (22) 

is superimposed at the upper left corner of each panel by 
a white curve. 

This means that physical condition in this region is not 
suitable for star-forming clouds and the ambipolar diffu¬ 
sion is not important in most cases. For some cases with 
high metallicity (^10' -^Ze) and low ionization parame¬ 
ter {t] = 0 or ().01), this boundary shifts toward smaller 
B into a region where the cloud VBx/uf[ = 1 is able to 
collapse. In such cases, the ambipolar diffusion can be 
important. 

For the Ohmic dissipation, whose rate depends only 
on the density, we can define dissipation density range 
where > 1- Fig. [3 shows the upper/lower bound 

of the dissipation density range by the Ohmic loss as a 
function of metallicity for different ionization parameters 
(indicated by lines with different colors). Except for the 
case of 7 ^ = 0, the gas becomes dissipative from higher 
density, i.e., the lower bound of the dissipation range in¬ 
creases, with decreasing metallicity for Z < O.IZq. This 
is because, for lower metallicity and smaller amount of 
the dust, the grains dominate the charge only at higher 
density (see Fig. |2|). On the contrary, without ionization 
source 77 = 0, the lower bound increases with metallic- 
ity for > IQ-'^Zq. This is because the charge is car¬ 
ried by the dust grains in this case even at low densities 


(nn lO^cm see Fig. |2| due to the rapid drop-off 
of the electron fraction. As a result, the larger amount 
of charge carriers is available for the larger metallicity 
and so the gas becomes less resistive. The same effect 
is also responsible for the upturn of the lower boundary 
around the solar metallicity in the cases with 77 > 0. On 
the other hand, the upper bound of the dissipation range 
remains almost the same for all models. Consequently, 
the dissipative range becomes narrower with decreasing 
metallicity. Eventually, in the case with 77 = 1 (77 = 10), 
for < 1 O“^Z 0 (< respectively) the dissipation 

range disappears. Such clouds do not experience the dis¬ 
sipative phase during the collapse. In addition, we can 
see that the dissipative phase is shorter in the density 
range for the cases with larger ionization parameter be¬ 
cause of the larger amount of charge carriers. 

4. SUMMARY AND DISCUSSION 

We have studied the coupling between the magnetic 
fields and the gas in the range of metallicity Z = 0 — IZq 
and ionization rate by cosmic rays/radioactive elements 
0 — 10 times the solar neighborhood value. Eor this pur¬ 
pose, we have calculated the thermal and chemical evo¬ 
lution of prestellar cloud cores with the Jeans scale by 
solving non-equilibrium chemical reactions for important 
charge carriers including charged grains. We have found 
that, for the Milky-Way like environment with metal¬ 
licity z© and the cosmic ionization rate 10 ^^s the 
magnetic flux at the scale of the core dissipates from the 
gas by the Ohmic loss, regardless of the field strength, in 
the density range 10^^cm“^ ^ '^h ^ 10^^cm“^. This dis- 
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Fig. 5. — Same as Fig[3l but for the ionization parameter rj = 1 (model 3). 


sipation is due to the enhanced resistivity in this range 
by the dominance of grains as charge carriers. With de¬ 
creasing metallicity and increasing ionization rate, more 
charge is carried by electrons and ions relative to grains 
and, as a result, the dissipation density range by the 
Ohmic loss becomes narrower. For metallicity less than 
10 “^ — 1 O“^Z 0 , depending on the ionization rate, the 
magnetic flux is always frozen to the gas at the Jeans 
scale. Hence, the magnetic field present in the prestel- 
lar core with such low metallicity will be taken into the 
forming protostar without dissipation. 

In what follows, we discuss implications of the 
magnetic dissipation for the dynamics of star-forming 
clouds, as well as for the small-scale dynamo action in 
the metal-free case. 

4.1. Implications for the dynamics of prestellar cores 

In the Milky-Way like environment, temperature in 
the prestellar core increases nearly adiabatically in the 
range 10 ^^cm“^ ^ ^ (see Fig. [Th). In 

this period, increasing pressure stops the dynamical col¬ 
lapse and a hydrostatic core, which is called the first core, 
forms at the center. The first core subsequently contracts 
only quasi-statically with accretion of matter from the 
envelope. Because of the slow radial velocity, the first 
core rotates many times before significant contraction. 
If the magnetic field is present and tightly coupled to 
the gas in this phase, the field lines are twisted enor¬ 
mously. The enhanced magnetic pressure/tension force 
launch the outflows from the first core, as well as ac¬ 


celerates the cloud collapse by transffering angular mo¬ 
mentum efficiently by the magnetic breaking. On the 
other hand, if magnetic fields dissipate before the first- 
core phase, neith er outflows nor the m agnetic breaking 
will be operative (jMachida et al.ll^QQSf ). 

In Fig. [3 the density ranges of the first core are indi¬ 
cated along with the Ohmic dissipation range for com¬ 
parison. The shaded region represents the range of adia¬ 
batic phase 7 > 7 ad, where the first core is formed. Note 
that the density range of the first core is independent 
of the ionization models. In the case of no ionization 
source (7 = 0 ), the dissipation zone extends in the den¬ 
sity range much wider than that of the first core for all 
metallicities. Hence, MHD effects such as the outflow 
launching or the magnetic breaking would not work in 
the first-core phase. Note that the first-core phase dis¬ 
appears at lowest metallicities, < IO^^Zq. The upper 
density bound of the first-core phase marks the onset of 
another dynamical collapse, which is induced by the ef¬ 
fective cooling by the H 2 dissociation. Eventually, for 
density ur ^ 10 ^^cm“^, following the completion of H 2 
dissociation, the temperature begins increasing adiabat¬ 
ically again. The collapse slows down and another hy¬ 
drostatic core, called the second core or protostar, forms. 
In this phase, the magnetic coupling is tight and so the 
MHD effects can be important, as long as some mag¬ 
netic fields are still present, for all metallicities and ion¬ 
ization rates. With 7 = 0.01 and Z <, O.lZr:^, i.e., i n 
the first-galaxy-like environment (e.g. IWise et al.ir2012f ). 
the magnetic dissipation still begins before the first core 
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Model 4 (ri=10) 
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Fig. 6 .— Same as Fig|3| but for the ionization parameter 77 = 10 (model 4). 
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Fig. 7.— Density range where the magnetic fields dissipate by 
the Ohmic loss for the gas with different metallicities. Two curves 
of the same color indicate the lower/upper bound of the range. 
Four colors correspond to different ionization parameters 77 = 0, 
0.01, 1, and 10. Vertical dashed lines represent the lower bound 
of the metallicities below which no dissipation regions are found. 
Also shown for comparison is the density range for the first-core 
phase (shaded region), where the protostellar outfiow is expected 
to be launched if significant magnetic fields are present. Note that 
the first-core density range does not depend on the ionization pa¬ 
rameter. 


formation. This means that the MHD effect in such an 
environment would be weaker than in the present-day 
cases. On the other hand, with ionization rate similar to 


the Milky-Way value {r] = 1), the magnetic dissipation 
begins almost at the same density as the onset of the 
first core phase Z < 0.1 Zq or even later than that for 
Z > O.IZq. In the latter case, the magnetic coupling is 
still tight in the early first-core phase, and thus rotation 
of the first core twists the field lines at its edge. Hence, 
above O.IZq, the magnetic effects are important on the 
dynamics in the first-core phase as is known from studies 
on the present-day star formation. Finally, with ioniza¬ 
tion rate higher than the local value {rj = 10, starburst 
galaxy-like environment), the dissipation zone becomes 
even narrower. The dissipation begins always after the 
first-core formation at any metallicity, and the dissipa¬ 
tion zone eventually disappears for Z < 10~^Zq. Hence 
the MHD effect is always important as long as the first 
core forms. 

We plan to carry out numerical MHD simulations of 
low-metallicity star formation by utilizing the resistivity 
obtained in this calculation as a future study. 

4.2. Implications for small-scale dynamo action in 
primordial gas 

Recently, amplification of magnetic fields by the small- 
scale dynamo during the first-star for mation has at¬ 
tract ed attention of some authors (e.g. iSchober et al.l 
l2Q12[) . This theory assumes that parental clouds of 
the first stars are highly turbulent, as in the present- 
day molecular clouds, down to the viscous scale much 
smaller than the Jeans length. In such a circumstance, 
the turbulent energy will be transferred to magnetic en¬ 
ergy within the eddy timescale of the turbulence, if the 
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magnetic field is tightly coupled to the g as (|Kazaiitsevl 
119681: iBrandenburg fc SubramanianI [ 200 ^ . A magnetic 
field at the smaller scale has shorter eddy timescale and 
is amplified faster than that at the larger scale. The 
energy of magnetic field at the viscous scale is thus am¬ 
plified most quickly. The amplified magnetic field at a 
smaller scale is supposed to inversely cascades to larger 
scales. Finally, this continues until the magnetic field 
reaches the equi-partition level even at the Jeans scale. 

Now, we examine the validity of the flux freezing in 
the primordial gas, which is the basic assumption of the 
small-scale dynamo theory. As shown in Fig. [3]i, for 
77 = 0 (i.e. the pristine environment) the magnetic field 
is tightly coupled to the primordial gas {Z = 0). This is 
true, however, only at the Jeans scale Rj. Here, let us 
estimate the drift velocity vbx at a smaller scale. This 
should be compared with the typical turbulent velocity 
'^turb at this scale. The drift velocity is inversely propor¬ 
tional to the length scale (eq. [T 8 |) : vbx oc at a scale 
aRj with a < 1. For the spectrum of the turbulent ve¬ 
locity, we assume the power-law form, TXturb oc R^. Here, 
the power index 'd equals to 1/2 for the Burgers tur¬ 
bulence and 1/3 for the Kolmogorov turbulence, respec¬ 
tively. From consideration that the turbulence is driven 
by the gravitational collapse motion at the Jeans scale, 
the typical turbulent velocity at this scale is roughly 
given by the free-fall velocity Tiff. Thus the turbulent 
velocity at the scale aRj is 

^turb — ^ 'Eff¬ 
using these relations, the ratio is 

'EBcc _ VBx{Rj) 

'E'turb OiRj Tiff 


(23) 

(24) 


and the dissipation condition at the scale aRj is given 
by 


VBxjRj) ^ ^^+1 
US 


(25) 


Note that the left hand side in the above inequality have 
been shown in Fig. [3]i. 

On the other hand, the magnetic field grows fastest 
at the smallest scale of the turbulence, the viscous scale 
i?vis defined by the relation 


^turb-^vis — 


(26) 


where z^vis denotes the kinematic viscosity. Combinig 
the relations (|23]) and (|26]) and using the Reynolds num¬ 
ber Re = TiffRj/z/vis 5 the viscous scale can be given as 


Rvis — J fte 


-l/i'd+l) 


i.e., a = Re 




for the vis¬ 


cous scale. Note that the Reynolds number is approx¬ 
imately Re ^ 1 Q^(T 2 H/Icm~^)y^ in the present calcu¬ 
lation (see also iSchober et al]l 2 Q 12 [) . Substituting the 
parameter a = at the viscous scale into the 

inequality (|25]) . we have the dissipation condition at the 
viscous scale: 


VBx{Rj) 

Tiff 


> Re“^ - 10“^ 


/ riB 
Vlcm“^/ 


(27) 


It is worth noting that this condition does not depend 
on the power spectrum of the turbulence. _ 

For instance, at ^ lcm“^ as assumed in ISchober et al.l 
(|2Q12f ). the magnetic field is tightly coupled to the gas 
for T;Bcc(Rj)/TTff < 10“^. As seen in Fig. [3] d, this cor¬ 
responds to B < 10“^G at ^ lcm“^. Therefore, the 
magnetic field can be amplified up to r\j 10 by the 
small-scale dynamo at this density. Note that this max¬ 
imal field strength is roughly equals to 0.1 Her, which 
is dynamically non-negligible. Below < 10 ^cm“^, the 
boundary of the dissipative region (eq. [27|) stays at 
B - 10-^ - 10-^G, since the inverse of the Reynolds 
number decreases as the density increases (eq[27|), can¬ 
celling the positive slope of the level curves (Fig. [3lJ, 
the contours just below Her). This means that, despite 
of the cloud contraction, the magnetic field in the core 
remains roughly at the constant level as a result of the 
dissipation. Eventually, at ^ 10 ^cm“^, the inverse of the 
Reynolds number becomes ^ 10“^^, coinciding with the 
level of the vertical contour and the dissipation condi¬ 
tion of eq. ([27|) is satisfied. Hence, the magnetic fiux 
at the viscous scale dissipates from the gas in a dense 
circumstance such as > 10 ^cm“^. 

In summary, the rapid amplification of the field at the 
viscous scale up to the level of 10 “^ — 10 “^ G seems to 
be plausible, if the turbulence assumed here is present 
in the cloud of < 10 ^cm“^.@ It is worth noting that the 
amplified level of the field strength could be dynamically 
important. 

A caveat about the small-scale dynamo amplification 
is that, the nature of turbulence at very small scales, in¬ 
cluding the inverse cascading, in cosmological minihalos 
has not yet been studied in detail. A future effort to 
tackle this problem would be rewarded. 


We appreciate fruitful discussions with Masahiro 
Machida, and the support by Ministry of Education, Sci¬ 
ence, Sports and Gulture, Grant-in-Aid for Scientific Re¬ 
search (G22540295:HS, B25287040:KO). 
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APPENDIX 


He re we present t he expressions for the chemical heating/cooling rate in the present calculations, which is the same 
as in lOmnk^ (|2QQQ[) . The cooling rate (per unit volume) associated with the H2 formation/dissociation can be written 
by using the related reaction rates as: 

Mchem,H2 =4.48eVn(H2) (n(H)/cdisi + n(R2)kdis2) 

( 4.2eV \ 

- 0.2eV + --^- - - n(H.)nukgr 

V l + HH/Hcr/ 

3.53eVn(H)n(e)/cH- + 1.83eVn(R)n(R^)k^+ + 4.48eV (n(H)^A:3bodyi + n(H2)n(H)^/c3body2) 

- - - ( 1 ) 

(1 + nn/iicr) 

where Ucr is the critical density for H2 deexcitation given by equation (23) of lOmnkail (|2QQQ[) . and we have used 
the reaction rate coefhcints of H2 formation on the grain surfaces(/cgr; reaction 23 in Omukai 2000), through 
process {kn-; 8), process {k^+; 10) the three-body reactions (febodyi and /c3body2; 19 and 20), and the collisional 

dissociation rates {kdisi and kdis 2 ] 13 and 21). In the above, n(X) denotes the number density of the species “X”, while 
flu is the H nuclei number density. The net cooling rate associated with the H ionization is 


P^chem,H+ 13.6eV 


^dn(H+) 


dt 


krecn{}l'^)n{e) — /ccR,Hn(H) - /ccR,H2^(H2)) 


( 2 ) 


where the second term comes from the assumption that the photons emitted in the radiative recombination does not 
contribute to the heating and the third and forth terms from that the CR ionization via 

H + CR^H++e (3) 

H2 + CR^H+ + H + e, (4) 

whose rate coefficients are A:cr,h and /ccr,H 2 5 respectively, does not contribute to the gas cooling. Note that we consider 
the heating associated with the CR ionization separately. Similarly, the cooling rates associated with the He and He+ 
ionization, respectively, are given by 


P^cliem,He+ — 24.6eV ^ 


^dn (He+) 
dt 

^dn(He++) 


+ ^rec,He+ n(He+)n(e) - fccR.Hen(He) ) . 
+ fcrec,He++’^(He++)n(e)^ , 


P'^cliem,He+”l 79.0eV 


dt 


(5) 














13 


where Kec,i{e+ is the rate coefficient for the radiative recombination of He+ (reaction 4 in Omukai 2000) and kcR,B.e 
is the rate coefficient for the He ionization by CR: 

He ffi CR He+ + e. (6) 

Summing up all the contributions above, we have the chemical cooling rate: 

P-^chem P'^chem,H2 P-^cherrijH “1“ P'^chem,H+ “1“ P'^chem,He'^'^ ’ (^) 


In our calculation, heating/cooling associated with the H2 formation/dissociation can be important, while the recom¬ 
bination/ionization contribution has little significance owing to the small ionization degree. 

For high enough density (nn > lO^^crn”^), we employ the net change of the number density of H+, He+, He++ and 
H2 (Anx) times the latent heat divided by At as the chemical cooling/heating rate per unit volume: 


, _dn(He+) _dn(He++) ..^^n(H2) 

pAchem = 13.6eV—+ 24.6eV—^ + 79.0eV—^ - 4.48eV- ^ ^ 


dt 


dt 


dt 


dt 


( 8 ) 






